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Abstract 

We investigate the effects of electron correlations on the ground state energy 
and the chemical potential of a droplet confined by a parabolic potential 
at high magnetic fields. We demonstrate the importance of correlations in 
estimating the transition field at which the first edge reconstruction of the 
maximum density droplet occurs in the spin polarized regime. 
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The electronic structure of a quantum dot has been of great interest both theoretically 
and experimentally. The chemical potential /i^, the energy required to add one electron 
to a dot, of a droplet of N electrons can be directly measured from the conductance peaks 
resulted from resonant tunneling of an electron through a quantum dot .011 The magnetic 
field induces spin and orbital angular momentum changes of the ground state, thus, the 
resulting B — N phase diagram has been a subject of many theoretical studies.iHI The rich 
dependence of \x n on the applied magnetic field provides information on electron-electron 
interactions, especially, in the spin polarized regime, because any transition in this regime 
results from many-body phenomena which cannot be explained by a single electron picture. 
In recent experiments, Klein and his coworkers! reported a divergent spin susceptibility, 



failed to predict the size of the magnetic field window in which the maximum density droplet 



In this paper, we examine the correlation effects in a droplet containig a small number 
of electrons and calculate the ground state energy and the chemical potential via an ex- 
act treatment of electron-electron interactions which includes Landau level mixings. Our 
calculations show that the region of magnetic fields in which the ground state of the spin 
polarized droplet is the MDD state is much reduced, compared with the HF calculations, 
demonstrating the importance of correlations. With the Landau level mixing, we also find 
that magnetizations are reduced and susceptibility peaks are shifted to higher magnetic 
fields. We find bumps in the chemical potentials with varying the magnetic field, which 
gives rise to a ripple in the /i — B diagram which was observed by single electron capacitance 
spectroscopy.0 

For simplicity, we first consider the Hamiltonian (H) for a spin polarized droplet of N 
electrons confined by a parabolic potential \m*tJ^r 1 . The generalization to an arbitrary spin 
case is straightforward. In the strong magnetic limit where all electrons are in the lowest 
Landau level, many-electron states in symmetric gauge with total angular momentum L can 
be written as! 



which was explained successfully 




theory 



(MDD) state! is the ground state. 
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1 *L,N=l[(z i -z j )P M (z 1 ,Z2,...,z N )e w^ lZkl \ (1) 

i<j 

where z^ — — iy^ gives the coordinate of the fcth electron and the effective magnetic 



length / is given by I 2 = 2 m*to(B) ' wnere u o(B) = y^o + ^c/^ an d to c is the cyclotron 
frequency. Pm(zi, ■ ■ ■ , z n) is an Mth order symmetric polynomial of N variables, where 
M[= L — N(N — 1)/2] is the excess angular momentum. For each angular momentum L, the 
energy can be obtained by diagonalizing the Hamiltonian matrix for all possible iV-particle 
basis functions ^l,n- Once the ground state | L, N) is obtained, the contribution from the 
Coulomb interaction to the ground state energy can be expressed, 

(L,N\J2—\L,N) = ^-C(L,N), (2) 

where e is the dielectric constant and C(L, N) is a dimensionless quantity which only depends 
on L and N. 

To study the electronic structure of the droplet, previous theoretical attempts were 
mostly relied on Hartree-Fock theoryB or exact calculations which use a basis set in the 
Hilbert space restricted to the lowest one or two Landau levelsB! In our calculations, instead 
of extending the basis states to those belonging to higher Landau level indices, which are 
computationally too expensive, we employ a symmetric Jastrow correlation factor,0 which 
will be multiplied to ^l,n so that electron correlations are effectively included through 
the mixing of higher Landau levels. We test various kinds of the Jastrow factor using a 
variational Monte Carlo method and find that a simple form of exp[— J2i<j u ( r ij)]i where 
u(r) = ar + br 2 , gives the lowest energy as shown in Table |I[ Throughout this work, how- 
ever, we set a = to reduce computational demands and to perform computations in an 
analytical fashion. 

With information on C(L,N) in Eq. (Q), the ground state energy in the variational 
many-body state, $l,tv(&) = e~ ^ l<jTt:i ^ l,n, can be evaluated in an analytical fashion. 
First, let us decompose the Hamiltonian into two parts, Hq and Vq, where 

Hn - B^J^ + ^n 2 ] - ^W, - E ^ (3) 



en — r-; 



(4) 



«<j - - i ■ i j i 

where the total spin angular momentum S z is N/2 for a spin polarized droplet. If the 



parameter b in $l,jv(6) is replaced by m* n °^ ^ , where Q = Ju^B) — NQ 2 , § L N (Q) will 
be the eigenstates of Then, the eigenvalues of Hq for the states of zero center-of-mass 

angular momentum (a general case will be discussed later) are 



Hn$L,N(ty = [E n (N) + Lh(Q - cu c /2)} $ LiJV (fi), 



(5) 



where 



1 m* 

En(N) = hu (B) + (N- l)fin - -g*—huj c S z . 

2 m,Q 



(6) 



The expectation value of the first term (an ad hoc harmonic potential) in Eq. @j can be 
evaluated analytically by using the raising (a J and and lowering (a~- and 6~) operators 
defined as,0 



a ij 
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m*Q ^,- - i{Pij,x ~ iPij, y ) 



m 



'Qqz-j - i(p ij>x + ip ij>y ) 



(7) 



(8) 



where a^ib^) is the hermitian conjugate of afj(bfj), Pij = Pi — pj = {pij^p^^^ and 



Zi — Zj. Then, it follows that 



(9) 



Considering the states of zero center-of-mass angular momentum, the polynomial Pm in Eq. 
(H) can be expressed only in terms of relative coordinates, and the variational many-body 
state is rewritten with the use of Zij = ( m ^ n ) l ^ 2 {. a tj + Kj)i 



®L,n(®) = II Z ij P M(J2 Z U, E Z 2i, • • • , E Z Ni)^v 
i<j i^l i^N 

= f-*-W J] a+P A/ (E 4, E 4, • • • , E <&)*t 

i<j i^l i^N 



(10) 
(11) 



where ty v denotes the vacuum state, i.e., a^ v = = 0, and is given by ^ 



Nm*u (B) -^2 _ m*n 



Y,i<j r ij] where R is the center-of-mass coordinate. Using the com- 



ex PL 2h xv 2Nh ^Kj ' ij 

mutation relations of afj and bfj, the expectation value of the ad hoc harmonic potential 



is 



i<j 



\ hQ 2 I 



2tt 
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[NL + N(N - 1)]. 



(12) 



For the Coulomb interaction energy in the variational state $l,n{Q), replacing I by I\ 



T, 



2m* ' 



in Eq. (0) leads to the expectation value, 

(<s> LiN (n)\j2—\$L, N m) = 

i<3 €r ^ 



-C(L,N). 



(13) 



From Eqs. (H), (|12|) , and (|i3|) , we finally get the total energy E L N for total angular 
momentum L and zero center-of-mass angular momentum by minimizing (^^(fi) | H \ 
®l,n(&)) over the variable Q; 

TiVt 2 e 2 
E L , N = minu{E a (N) + Lh(n - lu c /2) + —[NL + N(N - 1)] + — C(L, N)}. (14) 

z\i el n 

In general, the total angular momentum is given by L = M r + M c , where M r and M c 
represent the angular momentums of the relative and center-of-mass motions, respectively 
If M c 7^ 0, the total energy should include an additional energy, M c %[ujq(B) —u c /2], with L 
replaced by M r in Eq. (|H|). However, one can see that the states with M c ^ give higher 
energies for given N and B. Thus, the ground state energy of the iV-electron droplet at 
magnetic field B is a minimum value of El,n calculated for L > N ^ N ~^ an d M c = 0. 

For a droplet with arbitary spin states, we can generalize the total energy E L ^ N to E^s z ,n 
by replacing C(L, N) with C(L, S z , N), which is defined analogously to Eq. (|]). In this case, 
instead of *$?l,n in Eq- (0), we use ^l,s z ,n defined as, 



L,S Z ,N 



e 4i 



Zi, Z 2 , ■ ■ ■ , ZNjXni 



(15) 
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where is an Lth order polynomial of iV variables, Xn is a A^-product of single particle 
spinors, and n denotes possible spinor configurations. Although ^l,s x ,n is totally antisym- 
metric, each of and Xn does not necessary satisfy the antisymmetric property. 

The ground state energies in the variational many-body state $>l,s z ,n for a six-electron 
quantum dot are plotted as a function of magnetic field in Fig. |l| Throughout this work, we 
use the confinement strength of Tiujq = 2.5 meV and the usual values for e and m* in GaAs. 
Compared with the results of conventional exact calculations using ^l,s z ,n in Eq. (|T5]) , our 
calculated energies are generally lower with significant reductions at lower magnetic fields. 
This is because the Jastrow factor introduces more effectively the Landau level mixing and 
the correlation effect in the region of low magnetic fields, which are absent in the conventional 
scheme. From Fig. [I], one can notice that with the Jastrow factor magnetizations (—dE/dB) 
are relatively decreased whereas susceptibility peaks are shifted to higher magnetic fields. 
We also note that the variational ground states exhibit different orbital and spin angular 
momentum quantum numbers, compared with the conventional exact calculations. 

At low magnetic fields, < B < B c = 1.1 T, the quantum dot is found to favor the equal 
occupations of the up- and down-spin states, i.e., L = 6 and S z = 0, while the minimum total 
spin state only appears in a very narrow range near zero field in the conventional calculations 
without the Jastrow factor. At B c , the transfer of electrons between the up- and down-spin 
states begins and continues due to the Zeeman splitting energy up to the value of Bj = 
1.25 T, where a spin polarized state appears. With varying the magnetic field, the bumps 
in the chemical potentials (hn = - En_i), which are associated with the ground state 
transitions, are clearly seen in Fig. |2], exhibiting a ripple in the fijy — B diagram. In fact, 
such a feature in the hn — B graph was observed in single electron capacitance spectroscopy 
measurement s .0 

Our calculations show a spin polarized droplet (S z = N/2) above Bi, and a magnetic 
field window in which the ground state of the spin polarized droplet is the maximum den- 
sity droplet [L = N(N — l)/2] where electrons occupy one-electron eigenstates of angular 
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momentum indices, m = 0, 1, and (N — 1). As the magnetic field increases, since the 
Coulomb interaction energy grows due to the reduction of the radius of the MDD, the MDD 
must undergo a transition to a higher angular momentum state (L = 20 and S z = 2) to 
lower the ground state energy. In the N = 6 droplet, the edge of the MDD is reconstructed 
at B r = 3.7 T, and at this field an abrupt upward shift of the chemical potential appears, 
as shown in Fig. |[ Such a step at B r is the reminiscent of the chemical potential jump 
in the fractional quantum Hall effect, which is caused by quasi-hole creation^ Indeed, ex- 
perimentally, an upward step in /ijv at B r was observed for a droplet in the spin polarized 
regime.! 

Recently, Klein and his coworkers observed many transition fields in a droplet of N ~ 30 
electrons, and showed that the HF theory describes successfully their experimental data with 
a fitting parameter hu>o for B < However, the transition field B r was overestimated 
severely. Because of the computational difficulties for higher values of N, our calculations 
are restricted to droplets of N < 7. Nevertheless, we are able to see the effect of electron 
correlations on the transition fields. By varying the confinement strength from 1 to 5 meV 
for the N = 7 droplet, we plot the size of the magnetic field window, (B r — Bj), in which the 
MDD is the ground state of the spin polarized droplet, as a function of Bj in Fig. ^. With 
the Jastrow correlation factor, we find that the size of (B r — Bi) is significantly reduced by 
about 1 T, compared with the HF calculations, while in the conventional exact calculations 
with only the lowest Landau level the reduction is about 0.5 T. Thus, it is important to 
include the correlation effect in predicting the transition fields for the MDD state. 

In conclusion, we have made a detailed study of the correlation effect on the ground state, 
the chemical potential, and the transition fields in strong magnetic fields. Using an exact 
treatment of electron-electron interactions through the Landau level mixing, we find that 
magnetizations are reduced while the transition fields are shifted to higher magnetic fields 
by the correlation effect. We present a ripple in the /xjv — B diagram, which has not been 
explained by previous calculations. In estimating the size of the magnetic field window in 
which the MDD appears, we attribute the failure of the HF theory to the neglect of electron 
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correlations. 
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FIGURES 

FIG. 1. The ground state energies vs magnetic field for a six-electron quantum dot with the 
confinement strength of Tiujo = 2.5 meV. The arrows indicate the transition fields at which the 
ground state crossings occur. Some chosen ground states are labeled by the quantum numbers 
(L,2S g ). 

FIG. 2. (a) The chemical potential /ijy vs magnetic field for huJo = 2.5 meV. (b) The details 
of fiN for N = 7. Curve segments between two downward arrows are labeled by the ground 
state quantum numbers (L,2S Z ) of the N = 7 droplet, while those between two upward arrows 
correspond to the ground states of the N = 6 droplet. 

FIG. 3. For the spin polarized droplet of N = 7 electrons, (B r — Bj) is plotted as a function of 
Bj which is determined by varying the confinement strength from 1 to 5 meV. The present results 
with the Jastrow correlation factor are compared with the Hartree-Fock and conventional exact 
calculations without the Jastrow factor. 
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TABLES 

TABLE I. The variational ground state energies calculated by a Monte Carlo method are 
compared for various forms of the Jastrow factor, exp[— J2i<j u ( r ij)]- We use N = 3, B = 3 T, 
and Tiujq = 5.4 meV. 

u(r) Ground state energy (meV) 

br 2 41.041 
ar + br 2 40.866 ± 0.005 

alog(r) 41.248 ±0.010 

br/{l + ar) 41.230 ± 0.010 

^[1 -exp(- v /f/6-r/26)] 41.194 ± 0.005 

a[l - exp(-6r)]/r 41.750 ± 0.100 
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